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SECTION  H 

F-k  analyses  of  short-period  recordings  of  both  explosion  and  earthquake  sources  at  the  high-frequency  NORESS 
array  indicate  secondary  arrivals  fiom  a  near-rcccivcr  source  about  25-30  km  southwest  of  the  array.  Use  of  residual 
seismograms,  derived  by  subtracting  lire  beamed  record  from  each  array  channel,  improved  the  identification  of  the 
same  scattercr.  F-k  power  difference  plots,  obtained  by  subtracting  (after  normalization)  the  f-k  power  versus  slowness 
estimates  for  the  initial  P  window  from  those  for  the  later  windows,  provided  nearly  identical  results.  The  secondary 
seismic  source  appears  to  be  short-period  surface  waves,  Rg,  generated  by  the  scattering  of  incident  P  waves  in  the 
region  of  Lake  Mjosa,  27  km  SW  of  NORESS,  where  there  is  about  1  km  of  relief  from  the  bottom  of  the  lake  to  the 
top  of  an  adjacent  hill.  Polarization  analyses  of  two  explosions  recorded  on  three-component  elements  of  NORESS 
also  suggested  a  secondary  phase  with  polarization  characteristic  of  Rg  and  arrival  time  and  azimuth  consistent  with 
those  derived  from  the  f-k  analyses  for  one  but  not  for  the  other  shot.  Near-source  scattering  is  investigated  by  f-k 
analysis  of  data  from  reciprocal  arrays  formed  by  interchanging  the  roles  of  source  and  receiver.  Application  to 
closely  spaced  Yucca  Flat  (NTS)  shots  recorded  at  several  stations  provides  evidence  for  Rg  — ►  P  scattering  near  the 
northwestern  and  other  basin  boundaries  of  the  Yucca  Valley.  It  seems  that  deterministic  locations  of  near-receiver 
P  —>  Rg  and  near-source  Rg  -*  P  scattering  can  be  obtained  by  f-k  techniques. 


SUMMARY 


SECTION  A 

F-k  analyses  of  short-period  recordings  of  both  U.S.  and  Soviet  underground  nuclear  explosions 
at  the  Eskdalemuir  (EKA),  Scotland  array  indicate  arrivals  from  a  local  source  about  15  km  due  north¬ 
west  of  the  array.  Use  of  residual  seismograms  derived  by  subtracting  the  beamed  record  from  each 
array  channel  confirmed  the  presence  of  the  same  scattcrcr.  The  secondary  seismic  source  appears  to  be 
short-period  surface  waves  due  to  the  scattering  of  incident  P  waves  at  a  nearby  deep  valley  known  as 
Moffat  Water,  in  agreement  with  Key’s  (1967)  earlier  results.  Similar  analyses  of  the  NORESS  (Nor¬ 
way)  array  data  from  U.S.  and  Soviet  nuclear  shots  also  indicate  secondary  arrivals  suggesting  a  local 
near-surface  scatterer  about  25-30  km  south-west  of  the  array,  in  the  region  of  Lake  Mjosa.  Thus, 
Short-aperture  array  data  can  be  useful  in  identifying  and  locating  sources  of  near-receiver  scattering. 

SECTION  B 

F-k  analyses  of  short-period  recordings  of  both  explosion  and  earthquake  sources  at  the  high- 
frequency  NORESS  array  indicate  secondary  arrivals  from  a  ncar-rccciver  source  about  25-30  km 
southwest  of  the  array.  Use  of  residual  seismograms,  derived  by  subtracting  the  beamed  record  from 
each  array  channel,  improved  tire  identification  -of-tho- same  scatters!.  F-k  power  difference  plots, 
obtained  by  subtracting  (after  normalization)  the  f-k  power  versus  slowness  estimates  for  the  initial  P 
window  from  those  for  the  later  windows,  provided  nearly  identical  results.  The  secondary  seismic 
source  appears  to  be  short-period  surface  waves,  Rg,  generated  by  the  scattering  of  incident  P  waves  in 
the  region  of  Lake  Mjosa,  27  km  SW  of  NORESS,  where  there  is  about  1  km  of  relief  from  the  bot¬ 
tom  of  the  lake  to  the  top  of  an  adjacent  hill.  ^  Polarization  analyses  of  two  explosions  recorded  on 
three-component  elements  of  NORESS  also  suggested  a  secondary  phase  with  polarization  characteris¬ 
tic  of  Rg  and  arrival  time  and  azimuth  consistent  with  those  derived  from  the  f-k  analyses  for  one  but 
not  for  the  other  shot.  Near-source  scattering  is  investigated  by  f-k  analysis  of  data  from  reciprocal 


arrays  formed  by  interchanging  the  roles  of  source  and  receiver.  Application  to  closely  spaced  Yucca 
Rat  (NTS)  shots  recorded  at  several  stations  provides  evidence  for  Rg  — >  P  scattering  near  the 
northwestern  and  other  basin  boundaries  of  the  Yucca  Valley.  It  seems  that  deterministic  locations  of 
near-receiver  P  — »  Rg  and  near-source  Rg  — >  P  scattering  can  be  obtained  by  f-k  techniques. 
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BROADBAND  F  K  ANALYSIS  OP  ARRAY  DATA  TO  IDENTIFY  SOURChS  OF  LOC  AL  SCATTERING 
I.  N  Gupta,  ('  S  Lynnes.  And  R  A  Wagner 
Teledyne  Geotech.  Alexandria,  Virginia 


Abstract.  F-k  analyses  of  short-period  recordings  of  both 
U  S.  and  Soviet  underground  nuclear  exploeions  at  the  Esk- 
dalemuir  (EKA),  Scotland  array  indicate  arrivals  from  a 
local  source  about  15  km  due  northwest  of  the  array.  Use 
of  residual  seismograms  derived  by  subtracting  the  beamed 
record  from  each  array  channel  confirmed  the  presence  of 
the  same  scattercr.  The  secondary  seismic  source  appears  to 
be  short-period  surface  waves  due  to  the  scattering  of 
incident  P  waves  at  a  nearby  deep  valley  known  as  Moffat 
Water,  in  agreement  with  Key’s  (1967)  results.  Similar  ana¬ 
lyses  of  the  NORESS  (Norway)  anay  data  from  U.S.  and 
Soviet  nuclear  shots  also  indicate  secondary  arrivals  sug¬ 
gesting  a  local  near-surface  scatterer  about  25-30  km 
southwest  of  the  array,  in  the  region  of  Lake  Mjosa.  Thus, 
short-aperture  array  data  can  be  useful  in  identifying  and 
locating  sources  of  near-receiver  scattering. 

Introduction 

There  have  been  few  studies  of  teleseismic  P  coda  in 
which  actual  sources  of  scattering  have  been  identified 
Perhaps  the  most  convincing  study  so  far  has  been  the 
identification  of  scattered  arrivals  from  the  region  of  a  nar¬ 
row  lake  in  Scotland,  about  13  km  long  called  Moffat 
Water  which  lies  at  the  bottom  of  a  northeast  southwest 
trending  200  m  deep  valley  (Key,  1967).  By  beaming  Esk- 
dalemuir  (EKA)  array  records  from  several  earthquakes  and 
nuclear  explosions.  Key  identified  a  low-frequency.  Ray¬ 
leigh  wave  arrival  with  velocity  of  about  2.5  km/sec  from 
the  direction  315°,  originating  about  13  km  from  the  center 
seismometer.  Theoretical  studies  by  Hudson  and  Boore 
(1980)  confirmed  the  source  ot  the  local  wave  to  be  due  to 
scattering  of  incident  P  into  Rayleigh  waves.  In  this  study, 
we  demonstrate  the  usefulness  of  broadband  f-k  analysis  to 
identify  sources  of  Leal  scattering  such  as  Moffat  Water 
near  EKA  and  Lake  Mjosa  southwest  of  the  small-aperture 
NORESS  array.  Both  U.S.  and  Soviet  nuclear  explosions 
are  used  as  seismic  sources  to  provide  different  source- 
receiver  azimuths  Residual  seismograms  formed  by  „ub 
trading  the  array  beam  from  each  channel  improve 
identification  of  the  sources  of  locally  scattered  arrivals 

F-K  Analysis  of  U  S.  and  USSR.  Shots  at  fiKA 

Array  recordings  of  seismic  events  are  useful  for 
studying  detailed  characteristics  of  wave  propagation  across 
the  array  We  use  a  broadband  f-k  technique  described  by 
Kvaema  and  Ringdal  (1986) 


(7>pyri<)ht.  1  l*y  »}*•  i '  in  <  i 1 

nurtor  i ? 


The  method  was  first  applied  to  the  Eskdalemuir  (EKA) 
array  recordings  of  two  underground  nuclear  explosions. 
The  EKA  array  has  20  elements  arranged  in  an  "L" 
configuration  with  element  spacing  of  0.9  km  and  max¬ 
imum  element  spacing  of  9.8  km  (Bache  et  al.,  1985).  The 
results  are  expressed  in  terms  of  frequency-slowness  spec¬ 
tral  estimates  of  ihe  signal  at  specified  times.  Starting  a 
few  seconds  before  the  onset  of  P,  Parzen-tapercd  6.4  sec 
long  windows,  with  incremental  shifts  of  0.2  sec,  were  used 
to  obtain  the  frequency-slowness  plots.  A  single  element 
record  of  the  Nevada  Test  Site  (NTS)  shot,  KASSER1  (28 
October  1975,  m,,  =  6.2,  epicentral  distance,  A  =  72°, 
back-azimuth  309°)  is  shown  in  Figure  la.  The  results  of 
f-k  analysis  for  several  time  windows  are  shown  in  Figure 
2.  Here  the  P  window  is  defined  as  that  centered  on  the 
onset  of  P  so  that  it  contains  only  3.2  sec  of  signal  The 
window  starting  1  sec  after  the  P  window  (designated  the 
P+1  sec  window)  was  generally  found  to  have  Ihe  highest 
power  and  may  therefore  be  considered  to  represent  the  ini¬ 
tial  P.  On  each  slowness  plot,  the  highest  amplitude  value 
is  indicated  (with  power  in  db)  along  with  its  direction  of 
approach  or  back-azimuth  (in  degrees)  and  phase  velocity 
(in  km/sec).  The  contour  interval  is  I  db  and  a  range  of  10 
db  is  used.  The  signal/noise  ratio  for  the  EKA  recordings  of 
KASSER1  is  good  over  the  frequency  range  of  0,5  to  3.0 
Hz.  The  corresponding  broadband  results  for  the  P+1  sec 
and  P+5  sec  windows  are  shown  in  Figures  2a  and  2b, 
respectively.  Figure  2c  shows  results  for  the  narrow  fre¬ 
quency  passband  of  I ±0.25  Hz.  Dominant  energy  in  the 
first  window  shows  a  back-azimuth  of  305°  and  phase  velo¬ 
city  of  about  19  km/sec,  both  are  close  to  the  expected 
values  lor  initial  V  The  power  contours  arc.  as  expected, 
elongated  along  directions  that  coincide  with  the  (Elections 
of  the  two  arms  of  ihe  EKA  array.  In  addition  to  main 
arrivals  from  the  ei  plosion  source  region.  Figures  2b-c  sug¬ 
gest  low-phase-velocity,  secondary  sources  of  energy,  espe¬ 
cially  from  the  northwest  direction,  probably  due  to 
forward  scattering  of  teleseismic  P  to  surface  waves,  as 
suggested  b)  Key's  ('.'LV, 

In  order  lo  enhance  the  local  secondary  arrivals  with 
respect  to  the  primary  arrival,  residual  P  and  P  coda  were 
obtained  by  subtracting  the  lime-domain  beam  from  the 
original  P  waveform  for  each  EKA  channel.  The  beam  used 
for  this  purpose  was  obtained  by  aligning  the  first  peaks 
after  the  onset  of  P  and  summing  all  channels.  The  result¬ 
ing  residual  seismograms,  expected  lo  be  richer  in  near- 
rccetver  effects  since  the  principal  arrival  due  to  direct  P 
has  been  suppressed,  were  then  used  to  obtain  the 
frequency-slowness  plois  Results  for  the  P+5  sec  window 
but  for  three  different  frequency  passbands  (Figures  2d -0 
indicate  a  consistent,  low  phase  velocity  principal  arrival 
A  comparison  ol  Futures  2b  c  with  2d  c  demonstrates  the 
usefulness  ot  residual  seismograms  in  ihe  identification  and 
fixation  of  sources  of  near  receiver  scattering 
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Fig.  I.  Single  dement  records  of  (a)  KASSERI  at  EKA 
array  and  (b)  Kazakh  shot  of  16  December  1984  and  (c) 
LOCKNEY  at  NORESS  array.  Zero-to-peak  amplitudes  are 
given  in  nanometers. 


The  maximum  amplitude  arrivals  in  Figures  2d-f  indicate 
a  phase  velocity  of  about  3  km/sec  and  back-azimuth  of 
about  310°.  Combining  this  with  the  phase  velocity  and 
back-azimuth  in  Figure  2a  and  assuming  the  scattering 
takes  place  near  the  surface,  the  scatterer  should  tie  about 
13  km  from  the  center  element  of  the  array.  These  results 
agree  remarkably  well  with  those  of  Key  (1967).  For  a 
center  frequency  of  2  Hz  (Figure  2f).  the  distribution  of 


secondary  arrivals  appears  to  be  almost  random,  with  many 
sources  of  nearly  equal  amplitudes.  It  is  interesting  to  note 
that  Key  (1967)  found  the  spectrum  of  the  local  wave  to 
peak  around  I  Hz  and  drop  sharply  at  higher  frequencies. 

A  similar  analysis  was  applied  to  EKA  array  records  of 
Eastern  Kazakh  shot  of  23  December  1979  (m^  -  6.2,  A  = 
47°,  back-azimuth  61°).  The  most  prominent  secondary 
arrival  appears  in  the  P+7  sec  windows,  with  nearly  the 
same  direction  and  phase  velocity  as  that  for  the  U  S  shot 
This  arrival  time  is  later  than  those  in  Figures  2b-e,  which 
is  consistent  with  the  longer  travel  path  of  the  scattered 
incident  P  wave  relative  to  the  direct  P.  Assuming  near- 
surface  scattering,  these  results  indicate  the  scattering 
source  to  be  about  13  km  from  the  array  center,  in  the 
vicinity  of  Moffat  Water,  again  in  agreement  with  Key’s 
results. 

F-K  Analysis  of  U.S.  and  U.S.S.R.  Shots  at  NORESS 

We  now  show  results  from  similar  analysis  of  data  from 
the  short-aperture  array,  NORESS,  which  consists  of  sen¬ 
sors  within  a  3  km  diameter  array  (Ingate  et  al .  1985).  The 
24  short-pen od  vertical-component  sensors  used  in  this 
study  are  located  along  four  concentric  rings.  A  single  ele¬ 
ment  record  of  the  Eastern  Kazakh  shot  of  16  December 
1984  (ntfc  -  6.1,  A  =  38°,  back-azimuth  75°)  is  shown  in 
Figure  lb.  A  shorter  (3.2  sec  long)  window  is  used  for  f-k 
analysis  because  of  the  smaller  dimensions  of  NORESS  as 
compared  to  EKA.  Results  are  shown  in  Figure  3  in  which 
the  initial  P  window  indicates  the  largesl-ampUtude  arrival 


(d)  P  ♦  5  0  sec.  1  75  t  1  25  Hi 
At:  1111  VmBPw  12  ? 


(e)  P  ♦  5  C  sec.  1.00  ±  0  25  Hr  (I)  P  ♦  5  0  sec.  2  00  i  0  50  Hr 

M  313  9  V*  lit  P«  111 
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Fig.  2.  Frequency-slowness  spectral  estimates  derived  from  (a  c)  normal  and  (d  f)  residual  EKA  records  of  KASSliRI  The 
numbers  on  top  of  each  plot  indicate  back-azimuth  (deg),  phase  velocity  (km/sec),  and  signj  power  (dD),  respectively  of 
the  highest -amplitude  arrival  which  is  shown  connected  to  the  origin  The  maximum  slowness  value  is  0  4  scc/km  I  he 
signal  windows  and  frequency  passbands  are  indicated  on  each  plot  Note  the  low  phjse  velocity  arrival  (torn  the  northwest 
direction  appearing  as  a  secondary  arrival  in  plots  (b)  and  (c),  and  as  the  principal  arncal  in  plots  (d  f) 
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Fig  3  Similar  to  Fig.  2  but  derived  from  NORESS  records  of  the  Kazakh  shot  of  16  December  1984  A  low  phase  velo¬ 
city  arrival  from  the  southwest  direction  appears  as  the  principal  phase  on  all  five  later  arrival  windows. 
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The  low  phase  velocitv  aniv.il  from  the  southwest  direction  is  hardly  seen  in  fb).  appears  .is  .1  secondary  arrival  in  (c)  hut 
appears  as  the  principal  arrival  in  pirns  (d  ft 
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to  have  a  back  azimuth  of  80°  (Figure  3a),  in  good  agree 
ment  with  the  e  spec  ted  value.  Two  later  windows  with  fre¬ 
quency  passbands  of  4±3  Hz  and  210.5  (Figures  3b-c)  indi¬ 
cate  prominent  secondary  arrivals  from  the  southwest  dircc 
lion.  Figures  3d-f  show  the  residual  seismogram  f-k  results 
for  passbands  of  4 ±3  Hz,  1±0.25  Hz,  and  2±0.5  Hz.  Each 
of  these  figures  indicate  low -phase -velocity  arrivals  from 
the  southwest.  A  comparison  of  Figures  3d-e  with  3b< 
shows  that  the  use  of  residual  seismograms  again  facilitates 
the  isolation  of  secondary  sources.  Assuming  the  secondary 
source  to  be  near  the  surface,  the  results  in  Figure  3  sug¬ 
gest  the  scatterer  to  be  25-30  km  from  the  array.  It  is 
interesting  to  note  that  there  is  a  large  lake  (Lake  Mjosa), 
trending  northwest-southeast,  about  24  km  southwest  of 
NORESS.  There  is  steep  topography  surrounding  the  lake 
in  this  direction;  the  hills  immediately  to  the  west  of  the 
lake  have  an  elevation  of  600  m  above  the  lake,  which  is 
approximately  400  m  deep,  producing  total  relief  of  1000  m 
over  a  short  horizontal  distance.  It  seems  likely  that  the 
scatterer  is  associated  with  this  prominent  feature. 

Frequency-slowness  plots  derived  from  NORESS  records 
of  the  NTS  shot,  LOCKNEY  (24  September  1987,  m*  - 
5.7,  A  =  74°.  back-azimuth  319°)  are  shown  in  Figure  4. 
The  results  are  similar  to  those  in  Figure  3,  indicating 
nearly  the  same  azimuthal  direction  for  the  scatterer.  The 
scattered  arrivals  in  plots  derived  from  residual  seismo¬ 
grams  (Figures  4d-f)  arc  significantly  more  prominent  than 
those  from  normal  seismograms  (Figures  4b-c).  Again  the 
results  in  Figure  4  suggest  the  source  of  scattering  to  be 
about  25  30  km  from  the  array 

Conclusions 

F-k  spectral  analyses  of  short-period  data  from  (he  EKA 
and  NORESS  arrays  have  been  used  to  identify  and  locate 
sources  of  local  scattering.  Residual  seismograms  have  been 
found  to  be  useful  in  isolating  the  secondary  source  from 
the  effects  of  the  primary  source  Using  both  U  S.  and 
U  S  S  R  nuclear  explosions  as  (he  seismic  sources,  the 
scattered  arrivals  have  been  identified  as  low-velocity 


Kaylcigh  waves  originating  from  the  deep  valley,  Molfal 
Water  near  EKA  and  the  vicinity  of  Lake  Mjosa  near 

NORESS 
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SKCTION  1$ 
INTRODUCTION 


Few  studies  of  teleseismic  P  arrivals  have  been  successful  in  identifying  and  locating 
actual  sources  of  scattering.  One  of  the  most  convincing  analyses  has  been  the  identification 
of  scattered  arrivals  from  the  region  of  Moffat  Water,  a  narrow  lake  in  Scotland  about  8  miles 
long,  which  lies  at  the  bottom  of  a  600  ft  deep  valley  (Key,  1967,  1968).  The  valley  lies  in  a 
NE-SW  direction,  and  the  apparent  secondary  source  is  at  a  point  about  one  third  of  its  length 
from  the  northeastern  end,  where  it  cuts  through  a  basalt  and  dolerite  dyke.  By  beaming  Esk- 
dalemuir  (EKA)  array  records  from  several  earthquakes  and  nuclear  explosions.  Key  identified 
a  low-frequency  Rayleigh-wave  arrival  with  velocity  of  about  2.5  km/sec  from  the  direction 
31o°  and  originating  at  a  distance  of  about  13  km  from  the  center  seismometer  The  direction 
and  velocity  of  the  local  wave  did  not  change  in  spite  of  variation  in  source  azimuthal  direc¬ 
tion  of  over  50°,  suggesting  it  to  be  due  to  a  relatively  small  localized  heterogeneity. 
Particle-motion  analysis  of  three-component  data  confirmed  the  local  wave,  with  an  amplitude 
generally  about  20%  to  40%  of  that  of  the  incident  P  wave,  to  have  the  character  of  a  Ray¬ 
leigh  wave.  Theoretical  studies  by  Hudson  (1967)  and  Hudson  and  Boore  (1980)  also  indi¬ 
cated  the  source  of  the  local  wave  to  be  due  to  scattering  of  incident  P  into  Rayleigh  wave. 
Broadband  f-k  analysis  of  EKA  recordings  of  both  U.S.  and  Soviet  nuclear  explosions  further 
confirmed  the  secondary  seismic  source  to  be  from  Moffat  Water  (Gupta  et  al.,  1990).  Resi¬ 
dual  seismograms  obtained  by  subtracting  the  beamed  record  from  each  array  channel 
significantly  improved  identification  of  the  sources  of  locally  scattered  arrivals. 

Our  recent  f-k  analyses  of  data  from  the  high-frequency  seismic  array  NORESS  (Nor¬ 
way)  suggested  the  presence  of  a  local  scatterer  about  25-30  km  southwest  of  the  array,  in  the 
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region  of  Lake  Mjosa  (Gupta  et  al. ,  1989;  Gupta  et  al. ,  1990).  In  this  more  detailed  study, 
Soviet  and  U.S.  nuclear  explosions  and  an  earthquake  are  used  as  the  seismic  sources,  provid¬ 
ing  significantly  different  source-receiver  azimuthal  directions.  Two  new  f-k  techniques  based 
on  the  use  of  residual  seismograms  and  f-k  power  difference  plots  help  improve  identification 
of  the  scatterer.  Similar  results  are  obtained  by  polarization  analysis  of  three-component 
NORESS  data. 

In  simple  terms,  the  principle  of  reciprocity  means  that  the  source  and  receiver  positions 
in  a  seismic  experiment  can  be  exchanged  without  affecting  the  observed  seismograms.  Con¬ 
sequently,  single  observation  point  records  of  closely  spaced  seismic  sources  can  be  used  to 
form  an  array  on  which  conventional  array  processing  techniques  can  be  applied  (Spudich  and 
Bostwick,  1987).  Teleseismic  P  arrivals  from  Yucca  Flat  explosions  are  known  to  have  com¬ 
plex  waveforms  that  vary  considerably  from  one  shot  to  another.  Finite-difference  simulations 
have  explained  much  of  this  complexity  as  being  due  to  the  scattering  of  explosion-generated 
Rayleigh  waves  into  teleseismic  body  waves  because  of  laterally  varying  structure  (Mcl^iugh- 
lin  et  al. ,  1987;  Stead  and  Helmberger,  1988).  An  attempt  is  made  he  e  to  examine  near- 
source  scattering  for  Yucca  Flat  shots  by  f-k  analysis  of  data  from  several  stations  by  forming 
reciprocal  arrays.  The  results  are  interpreted  on  the  basis  of  the  known  geology  of  the  Yucca 
Flat  test  site.  Both  near-receiver  and  near-source  scattering  play  important  roles  in  the  moni¬ 
toring  of  underground  nuclear  tests.  For  example,  the  Lg  phase  from  explosions,  known  to  be 
strongly  influenced  by  scattering  (c.g.,  Gupta  and  Blandford,  1983),  is  often  used  for  the 
detection,  source  discrimination,  and  yield  determination  of  underground  nuclear  explosions. 
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F-K  AND  POLARIZATION  ANALYSIS  OF  USSR  AND  U.S.  SHOTS  AT  NORESS 


Array  recordings  of  seismic  events  are  useful  for  studying  detailed  characteristics  of 
wave  propagation  across  the  array.  Broadband  f-k  methods  have  been  used  to  detect  various 
seismic  arrivals  and  determine  their  phase  velocity  and  azimuth  (Kvaema  and  Doornbos, 
1986;  Kvaema,  1989;  Ringdal  and  Kvaerna,  1989).  The  f-k  technique  used  here  is  described 
by  Kvaema  and  Doornbos  (1986)  and  Kvaema  and  Ringdal  (1986). 

We  first  applied  the  broadband  f-k  spectra  method  to  the  NORESS  recordings  of  two 
underground  nuclear  explosions.  The  NORESS  small-aperture  array  consists  of  sensors  within 
a  3  km  diameter,  and  the  24  short-period  vertical-component  sensors  used  in  this  study  are 
located  along  four  concentric  rings  (Kvaerna,  1989).  The  results  are  expressed  in  terms  of 
frequency-slowness  spectral  estimates  of  the  signal  along  various  azimuthal  directions.  Start¬ 
ing  a  few  seconds  before  the  onset  of  P ,  Parzen-tapered  3.2  sec  (128  point)  windows,  with 
incremental  shift  of  0.2  sec,  were  used  to  obtain  the  frequency-slowness  plots.  Results  from 
the  Kazakh  shot  of  16  December  1984  (mb  =  6.1,  A  =  38°,  back-azimuth  75°)  are  shown  in 
Figure  1.  Here  the  P  window  is  centered  on  the  onset  of  P  so  that  it  contains  only  1.6  sec  of 
signal.  The  window  starting  1  sec  after  the  P  window  (designated  the  P  +  1  sec  window)  was 
generally  found  to  have  the  highest  power  and  may  therefore  be  considered  to  represent  the 
initial  P .  On  each  slowness  plot,  the  highest  amplitude  value  is  indicated  (with  power  in  db) 
along  with  its  direction  of  approach  or  back-azimuth  (in  degrees)  and  phase  velocity  (in 
km/sec).  The  contour  interval  is  1  db,  and  a  range  of  10  db  is  used.  The  signal/noisc  ratio  for 
the  NORESS  recordings  of  this  explosion  is  good  over  the  frequency  range  of  about  0.5  to  7.0 
Hz.  The  broadband  (1  to  7  Hz)  results  for  the  P  +  1  sec  and  P  +  11.4  sec  windows  are 


7 


(g)  (P  +  1 1 .4  sec)  -  (P  +  1 .0  sec)  (h)  (P  +  1 1 .2  sec)  -  (P  «•  1 .0  sec)  (i)  (P  +  1 1 .4  sec)  -  (P  +  1 .0  sec) 
Az:  234.0  Vel:  3.22  Pwr:  54.3  Az:  233.8  Val:  3  18  Pwr:  59.9  Az:  233  0  Vel:  3  08  Pwr:  59.0 


Figure  1.  Frcqucncy-slowncss  spectral  estimates  derived  from  NORFSS  records  of  the  Kazakh 
shot  of  16  December  1984  based  on  the  use  of  (a,  b,  c)  normal,  (d,  c,  0  residual  seismograms,  and 
(g,  n,  i)  f-k  power  difference  plots.  The  numbers  on  top  of  each  plot  indicate  hack-azimuth  (deg), 
phase  velocity  (km/scc),  and  signal  power  (dB),  respectively  of  lire  highcsl-ampiitudc  arrival  which 
is  shown  connected  to  the  plot  origin.  The  maximum  slowness  value  is  0.4  sec/km.  The  signal 
windows  and  frequency  passbands  are  indicated  on  each  plot.  Note  the  low  phase  velocity  arrival 
from  the  southwest  direction  appearing  as  the  principal  phase  on  all  eight  later  arrival  windows. 


shown  in  Figures  la  and  lb,  respectively.  Figure  1c  shows  results  for  the  narrow  frequency 
passband  of  2.0  ±  0.5  Hz.  Dominant  energy  in  the  first  window  shows  a  back-azimuth  of 
about  80°  and  phase  velocity  of  about  16  km/sec;  both  are  close  to  the  expected  values  for 
initial  P.  In  addition  to  main  arrivals  from  the  explosion  source  region,  Figures  lb  and  lc 
suggest  low-phase-velocity  secondary  sources  of  energy,  especially  from  the  southwest  direc¬ 
tion. 

The  local  secondary  arrivals  were  enhanced  with  respect  to  the  primary  arrival  by  sub¬ 
tracting  the  time-domain  beam  from  the  original  P  waveform  for  each  array  channel  and 
using  the  resulting  residual  seismograms  for  f-k  analysis  (Gupta  et  al. ,  1990).  This  simple  but 
effective  method  constitutes  a  first  step  in  iterative  beamforming  of  two  simultaneous  signals 
and  provides  a  first  approximation  to  a  maximum  likelihood  estimate  (Blandford  et  al. ,  1976). 
Results  for  the  same  two  windows  and  frequency  passbands  as  in  Figures  lb  and  lc  are 
shown  in  Figures  Id  and  le,  respectively;  a  comparison  clearly  demonstrates  the  usefulness  of 
residual  seismograms  in  the  identification  and  location  of  sources  of  near-receiver  scattering. 
Figure  If  shows  results,  based  on  the  use  of  residual  seismograms,  for  the  lower  frequency 
passband  of  1 .00  ±  0.25  Hz.  The  three  Figures  Id,  le,  and  If  consistently  indicate  a  prom¬ 
inent  low  phase  velocity  arrival  from  nearly  the  same  azimuthal  direction. 

An  attempt  was  also  made  to  investigate  secondary  sources  of  scattering  with  the  help  of 
f-k  power  difference  plots  obtained  by  subtracting  the  normalized  two-dimensional  matrix  of 
f-k  power  versus  slowness  estimates  of  the  initial  P  window  from  those  for  the  later  win¬ 
dows.  Results  for  the  same  three  frequency  passbands  as  used  in  Figures  Id,  le,  and  If  are 
shown  in  Figures  lg,  lh,  and  li,  respectively.  The  three  figures  indicate  the  same  prominent 
low  phase  velocity  arrival  from  nearly  the  same  azimuth.  It  appears  therefore  that  f-k  power 
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difference  plots  provide  another  simple  and  effective  method  for  studying  the  characteristics 
of  secondary  arrivals. 

* 

The  maximum  amplitude  arrivals  in  all  eight  figures  based  on  later  time  windows  (Fig¬ 
ures  lb  through  li)  indicate  a  phase  velocity  of  about  3  km/sec  and  back-azimuth  of  about 
235°  .  Combining  this  with  the  phase  velocity  and  back-azimuth  in  Figure  la  and  assuming 
the  scattering  to  take  place  near  the  surface,  the  scatterer  should  lie  25-30  km  from  the  array. 

There  is  a  large  lake,  Lake  Mjosa,  running  along  the  northwest  to  southeast  direction,  about 
24  km  southwest  of  NORESS.  The  hills  immediately  to  the  west  of  the  lake  have  an  elevation 
of  600  m  above  the  lake,  which  is  approximately  400  m  deep,  meaning  that  there  is  a  varia- 
don  in  elevation  from  the  bottom  of  the  lake  to  the  top  of  the  hill  of  1000  m  over  a  short  hor¬ 
izontal  distance  (S.  Mykkeltveit,  personal  communication,  1989).  It  seems  likely  that  the 
scatterer  is  associated  with  this  prominent  geological  feature  which  has  the  largest  topographi¬ 
cal  relief  within  at  least  50  km  of  the  array. 

Frequency-slowness  plots  derived  from  NORESS  records  of  the  NTS  shot  LOCKNEY 
(24  September  1987,  mb  =  5.7,  A  =  74°,  back-azimuth  319°)  are  shown  in  Figure  2,  which 
includes  results  from  normal  f-k  analysis  (Figures  2a,  2b,  and  2c),  from  residual  seismograms 
(Figures  2d,  2e,  and  20,  and  from  f-k  power  difference  plots  (Figures  2g,  2h,  and  2i).  Results 
are  similar  to  those  in  Figure  1,  indicating  nearly  the  same  azimuthal  direction  for  the 
scatterer.  The  secondary  arrival  is  most  prominent  in  the  P  +  10.6  sec  window.  The  earlier 
arrival  time,  as  compared  to  those  in  Figures  lb  through  li,  is  consistent  with  the  somewhat  » 

shorter  travel  path  of  the  scattered  incident  P  wave  relative  to  the  direct  P.  The  scattered 
arrivals  in  plots  derived  from  residual  seismograms  (Figures  2d,  2e,  and  20  as  well  as  those 
from  difference  f-k  plots  (Figures  2g,  2h,  and  2i)  are  significantly  more  prominent  than  those 
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(d)  P  +  10.6  sec.  1.25  ±  0.75  Hz  (e)  P  +  10.6  sec.  2.00  ±  0.50  Hz  (f)  P  +  10.6  sec,  1.00  ±  0.25  Hz 

Az:  244.9  Vel:  2.66  Pwr3.0  Az:  239.5  Vel:  2.86  Pwr  6.5  Az:  258.3  Vel:  3.00  Pwr:  0.3 


(g)  (P  +  10.6  sec)  -  (P  +  1.0  sec)  (h)  (P  +  10.6  sec)  -  (P  +  1.0  sec)  (i)  (P  +  10.6  sec)  -  (P  +  1.0  sec) 

Az:  268.6  Vel:  4.61  Pwr:  14S  Az:  248.5  Vel:  3.13  Pwr  23  8  Az:  253.4  Vel:  2.40  Pwr:  6  6 


Figure  2.  Similar  to  Figure  1  but  derived  from  NORESS  records  of  LOCKNEY  with  a  maximum 
slowness  of  0.5  sec/km.  The  low  phase  velocity  arrival  from  the  southwest  direction  is  hardly  seen 
in  (b),  appears  as  a  secondary  arrival  in  (c)  but  appears  as  the  principal  phase  in  nearly  all  other 
later  arrival  windows. 


from  normal  seismograms  (Figures  2b  and  2c).  Again  the  results  in  Figure  2  suggest  the 
source  of  scattering  to  be  about  25-30  km  from  the  array. 

Polarization  analyses  of  the  four  3-component  short-period  elements  of  NORESS  for  the 
USSR  and  U.S.  explosions  described  earlier  were  carried  out  using  the  method  of  Jurkevics 
(1988).  A  zero-phase  Butterworth  bandpass  filter  is  applied  to  the  data,  and  covariance 
matrices  are  averaged  over  the  array  for  a  sliding  window.  The  polarization  ellipsoid, 
representing  the  best  fit  to  the  data  in  a  least- squares  sense,  is  computed  from  the  average 

covariance  matrix.  The  characteristics  of  ground  motion  can  be  specified  in  terms  of  attri¬ 

butes  derived  from  the  principal  axes  of  the  polarization  ellipsoid.  Denoting  the  latter  by  Xj , 
hi,  and  X.3,  where  Xj  >  hi  ^  X,3,  the  degree  of  rectilinearity  is  given  by  ( 1 — (X2+X3 )/2X j )  which 
is  1.0  when  there  is  only  one  nonzero  eigenvalue,  as  for  pure  P  waves.  Pure  Rayleigh- wave 
motion  is  elliptical  and  the  panicle  motion  is  confined  to  a  plane.  Defining  oblateness  as 
2(X7-X3)(A.i-A.3)/(X.12+^2+X32),  Rayleigh  waves  should  be  characterized  by  very  oblate  (i.e.,  fiat 
elliptical)  polarization  (the  maximum  oblateness  is  1 ,  when  ^  and  X*  =  0).  Results  from 
the  Kazakh  shot,  based  on  the  0.5- 1.0  Hz  band,  are  shown  in  Figure  3,  indicating  a  prominent 
phase  with  a  polarization  characteristic  of  Rayleigh  waves  arriving  about  1 1  sec  after  P.  The 

back-azimuth  computed  from  the  polarization  is  about  220-250°.  The  arrival  time  and 

azimuth  are  therefore  consistent  with  those  for  the  secondary  seismic  source  in  the  f-k  ana¬ 
lyses  (Figure  1).  Similar  analysis  was  carried  out  for  the  NTS  shot  LOCKNEY  (Figure  4). 
The  rectilinearity  is  generally  high,  indicating  predominantly  P- wave  motion,  probably 
representing  telescismic  P  waves  from  the  source  region.  This  is  supported  by  the  f-k  ana¬ 
lyses  in  Figures  2b  and  2c,  which  indicate  a  large  amount  of  energy  arriving  form  the  source 
region  several  seconds  after  the  direct  P.  Rayleigh  wave  motion  is  not  apparent  <n  the 


12 


NHAO 


nhc:> 


'W'V - 'N/VV'-- 


W\Aaa^ 


r - t"  i  r - 1 - 1 - f  i  -» 

0  0  2.5  5  0  7.5  10  0  12  5  15  0  17.5  20  0 

NRC4  TIME  (sec) 


-  ^WV' — Wa--- 


AA/VW-- 


0.0  2.5  5  0  rs  10.0  17.5  15.0  17.5  20.0 


TIME  (sec) 


+  +  +  +  ♦  *  *  *  (e) 


^  0  •" 
E 

2i  °  s~ 

z 


MIW 


I - |— 7  — ■ — I - 1 - 1 - 1 - 1 - 1 

0  0  2.5  5  0  7  5  10.0  12.5  IS  O  17.5  20.0 


NRC7  TIME  (sec) 

.  A  .  (d) 


a/\W\Aa^ 


i - r - r — — i - 1 - 1 - 1  r  t 

0.0  2.5  5.0  7.5  10.0  12.5  IS  O  1  7.5  20.0 


TIME  (sec) 


o  © 


O 

UJ 

CC  o  2- 


©  O  Q  ©  ® 

o .  o~t - t - r - 1 - t - 1  t —  1  r  o.oi - r  i 


0.0  2.5  5  0  7.5  10  0  12  5  15.0  1  7.5  20  0 


TIME  (sec) 


0.0  2.5  5  0  7  5  10.0  12.5  15.0  1  7.5  20.0 

TIME  (sec) 


»ooooo»%(g) 


CO  o 

UJ 

Z3 

— !  0  6  - 
< 

>  * 

z 

UJ  0.4-  + 

C3 

UJ 

o  ?-  J 


O  o  ° 
o  o 

O  o  100- 

**  4*  F  5°- 

.  **  i  ° 


100  o  o  o  o  o  o 


0  0)-  T  T  I  -  T  T  I  1  1  -  1 

00  7  5  50  7  5  1  0  0  1  ?  S  1  S  0  1  7  5  ?0  0 


|  1  “  1  1  1 - 1 - 1 - 1 - 1 

00  75  50  75100  17  5  1 5  0  1 7  5  70  C 


TIME  (SCC) 


TIME  (sec) 


Figure  3.  Polarization  analysis  of  16  Dec  19X4  Shagan  River  event,  (a,  h,  c,  d)  Three-component 
NORESS  seismograms,  (c)  Recti  linearity  versus  time.  (0  Oblateness  versus  time,  (g)  Eigenvalue 
amplitudes  (normalized  to  l).  (ti)  Back-azimuth,  assuming  Rayleigh  wave  polarization:  line  shows 
back-azimuth  to  Shagan  River. 
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Figure  4.  Similar  to  Figure  3  for  the  NTS  shot,  LOCKNFY. 


polarization  until  about  14  seconds  after  the  direct  P.  The  back-azimuth  is  260-280°,  in  rough 
agreement  with  the  f-k  results. 
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F-K  ANALYSIS  OF  AN  INTERMEDIATE-FOCUS  EARTHQUAKE  AT  NORESS 


We  carried  out  f-k  analysis  on  NORESS  data  from  an  earthquake  in  Romania  (45.82°  N, 
26.65°  E,  August  1,  1985,  mb  =  4.7,  depth  122  km,  A  =  17°,  back-azimuth  142°)  in  order  to 
examine  any  possible  influence  of  variations  in  source-type  and  source-receiver  azimuth  on 
the  determination  of  near-receiver  scattering.  The  results  are  shown  in  Figure  5.  The  P  +  1 
sec  window  (Figure  5a)  indicates  the  expected  phase  velocity  and  back-azimuth,  whereas  the 
two  P  +  11.6  sec  windows  (Figures  5b,  c)  contain  some  evidence  of  a  low-phase- velocity 
secondary  arrival  from  the  southwest  direction.  Comparing  these  with  the  results  from  explo¬ 
sions  (Figures  lb,  c  and  2b,  c),  the  earthquake  P  coda  has  considerably  more  energy  coming 
from  the  source  region,  probably  because  of  the  significantly  longer  source  duration. 
Nevertheless,  f-k  based  on  residual  seismograms  (Figures  5d,  e,  f)  and  f-k  power  difference 
plots  (Figures  5g,  h,  i)  produce  results  for  the  secondary  source  that  are  remarkably  similar  to 
those  from  the  two  explosions  (Figures  1  and  2). 


16 


(a)  P  +  1.0  sec.  4.00  t  3.00  Hz  (b)  P  ♦  1 1  6  sec,  4.00  ±  3  00  Hz  (c)  P  +  1 1.6  sec,  2.00  +  0.50  Hz 

Az  145  9  Vet  10  27  Pwr  105  Az:1402  Vot:  I0.G4  Pwr:  10.7  Az:  138.3  Vet  9.82  Pwr:  14.9 


(d)  P  +  1 1 .6  sec.  4.00  ±  3.00  Hz  (e)  P  4-  1 1 .6  sec.  2.00  ±  0.50  Hz  (f)  P  +  1 1 .6  sec.  1 .00  ±  0.25  Hz 

Az:  227.0  Vet:  2.91  Pwr:  6.4  Az:  224.2  Vet:  3.06  Pwr:  9.9  Az.  228.6  Vet:  2.83  Pwr  8  5 


(g)  (P  +  1 1  6  sec)  -  (P  4  1.0  sec)  (h)  (P  +  1 1 .6  sec)  -  (P  +  1.0  sec)  (i)  (P  +  1 1 .6  sec)  -  (P  +  1.0  sec) 

Az:  229  6  Vet  2  81  Pwr  2  0  Az  225  1  Vet:3l7  Pwr  7  6  Az  227  7  Vet  289  Pwr  192 


l  tgurc  5.  Similar  to  ligure  1  but  derived  from  NORIiSS  records  of  die  Romanian  earthquake  of  1 
August  1085.  The  low  phase  velocity  arrival  from  the  southwest  direction  appears  as  a  secondary 
arrival  in  (b)  and  (c)  but  appears  as  the  principal  phase  in  all  other  later  arrival  windows. 
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F-K  ANALYSIS  OF  DATA  FROM  RECIPROCAL  ARRAYS 


We  also  used  f-k  analysis  to  investigate  near-source  scattering  by  forming  reciprocal 
arrays  in  which  the  roles  of  source  and  receiver  are  interchanged.  Short-period  records  at  a 
single  station  of  seismic  sources  that  are  fairly  similar  and  closely  spaced  may  be  used  for  this 
purpose.  In  order  to  simulate  a  reciprocal  array,  arrival  times  of  the  first  peak  in  the  direct  P 
wavetrain  from  each  explosion  are  used  as  reference  points,  and  the  traces  are  shifted  on  the 
basis  of  their  predicted  travel  times  to  the  recording  station.  Yucca  Flat  (NTS)  explosions 
recorded  at  several  stations  at  teleseismic  distances  were  used  to  construct  four  reciprocal 
arrays,  and  the  sources  of  secondary  arrivals  were  examined.  The  locations  of  39  shots  used  in 
each  of  the  four  reciprocal  arrays  used  in  this  study  are  shown  in  Figure  6.  Two  of  these 
arrays  were  derived  from  groups  of  10  shots  each  recorded  at  the  ASRO  station  MAJO,  these 
have  been  designated  as  MAJO-North  and  MAJO-South.  The  center  element  NAO  records  of 
the  NORSAR  array  for  15  shots  provided  another  reciprocal  array.  Lastly,  18  records  from  a 
combination  of  the  WWSSN  stations  TOL,  PTO,  and  MAI.  having  nearly  the  same  distance 
and  back-azimuth  to  Yucca  Flat  were  combined  to  form  the  fourth  reciprocal  array.  The  long 
dimensions  of  these  four  arrays  vary  from  about  3  to  5  km.  The  seismic  velocities  in  the 
Yucca  Flat  region  are  known  to  be  considerably  smaller  than  those  in  the  shield  region  of 
NORESS.  For  this  reason,  6.4  sec  windows  were  used  for  the  f-k  analysis,  and  the  maximum 
slowness  was  increased  to  1.0  scc/km  on  all  slowness  plots. 

Results  from  MAJO-North  array  are  shown  in  Figure  7.  Frequency-slowness  plots, 
obtained  with  an  incremental  shift  of  1  sec  for  three  different  frequency  ranges,  showed  a 
prominent  secondary  arrival  in  the  P  +  4  sec  window.  Frequency  slowness  plots  of  the  initial 
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TOL  (8)  81° 
PTO  (5)  78° 


Figure  6.  Locations  of  39  Yucca  Flat  explosions  used  in  forming  four  reciprocal  arrays.  (1) 
MAJO-Norlh  consists  of  10  shots  denoted  by  "J",  (2)  MAJO-South  with  10  shots  .silhin  the  en¬ 
closed  area,  (3)  NAG  consists  of  15  shots  denoted  by  "o",  and  (4)  TOL-PTO-MAL  3rray  consist¬ 
ing  of  15  records  of  1 1  shots.  Arrows  indicate  azimuthal  direction  to  each  station  or  mean  azimuth 
to  group  of  stations.  The  first  number  for  each  array  (in  parentheses)  indicates  die  number  of 
shots  and  the  second  the  distance  in  degrees. 
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(a)  P  +  10  sec.  1  25  i  0.75  Hz  (b)  l>  .  4  0  sec.  1  25  i  0  /5  Hz 

Az:  311.4  Vcl  1909  Pwr  14  3  A/  2/4  1  Vd:  5.73  Pm  19  3 


(c)  P  +  1.0  sec,  1.00  ±  0.50  Hz  (d)  P  +  4.0  sec,  1.00  ±  0.50  Hz 

Az:  309  5  Vet:  19.16  Pwr  : -14.1  Az:  275  7  Vel:  5  87  Pwr: -18  8 


(e)  P  +  1.0  sec,  0.50  ±  0.30  Hz  (f)  P  +  4  0  sec.  0.50  ±  0.30  Hz 

Az:  307.8  Vel:  25.52  Pwr: -14.5  Az  297  3  Vrt  8  37  Pwr: -20  4 


11)  1  BO  0  10  10  1  MO  il  10 


Figure  7.  Similar  to  Figure  1  bul  derived  from  the  reciprocal  array  MAJO-North.  The  maximum 
slowness  value  is  1.0  scc/km  and  results  for  three  dillerent  frequency  passhands  arc  shown.  The 
prominent  secondary  arrival  in  the  P  +  4  sec  windows  has  a  phase  velocity  of  about  1.6  km/sec 
and  azimuth  of  about  3(X)°. 
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P  (i.c .,  the  P  +  1  sec)  window  indicate  that  the  dominant  energy  arrives  close  to  the  expected 
azimuth  and  phase  velocity.  However,  unlike  the  results  from  the  initial  P  window  from 
NORHSS  (Figures  la,  2a,  5a),  there  are  several  other  secondary  sources  of  energy.  This  con¬ 
tamination  is  at  least  in  part  due  to  the  irregular  distribution  of  shot  locations  and  variations  in 
their  source  functions.  The  prominent  secondary  arrival  has  a  phase  velocity  of  about  1.6 
km/sec  and  arrives  with  an  azimuth  of  about  300°. 

Figure  8  shows  results  from  MAJO-South  array  for  the  same  three  frequency  ranges  as 
in  Figure  7.  A  prominent  secondary  arrival  appears  in  the  P  +6  sec  windows,  with  phase 
velocity  of  about  1.6  and  azimuth  of  about  260°.  Results  from  the  NAO  array  are  shown  in 
Figure  9.  The  secondary  arrival  appears  best  in  the  P  +7  sec  window  and  has  a  phase  velo¬ 
city  of  about  1.8  km/sec  and  azimuth  of  about  270°.  The  array  made  by  combining  18  records 
of  11  shots  at  3  stations  provided  results  shown  in  Figure  10.  Here  the  P  +7  sec  window 
shows  a  low  phase  velocity  (about  2.3  km/sec)  arrival  as  the  primary  phase,  with  azimuth  of 
about  117°.  This  arrival  in  particular  appears  very  similar  in  location  and  delay  time  to  an 
arrival  obtained  by  Lynnes  and  Lay  (1989)  from  WWSSN  and  Canadian  data  using  a  sem¬ 
blance  technique. 
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(a)  P  +  1.0  sec,  1.25  t  0.75  Hz 

Az:  282.9  Vd.  16  04  Pwr  14  6 


(b)  P  4  6  0  sec.  1.25  i  0.75  Hz 

Az  24/ ‘j  Vd  6.35  Pwr: -23.1 


90  0 


(c)  P  +  1 .0  sec,  1 .00  ±  0.50  Hz 

Az:  233.0  Vel:15.66  Pwr:-142 


9  6 

(d)  P  +  6.0  sec,  1.00  ±  0.50  Hz 

Az:  244, 1  Vd  6.33  Pwr:  -22.6 


10 


90.0 


-1.0 


(e)  P  +  1 .0  sec.  0.50  ±  0.30  Hz 

Az:  271.7  Vd.  12.15  Pwr: -14.7 


(1)  P  +  6.0  sec,  0.50  ±  0.30  Hz 

Az:  226  9  Vd  9  34  Pwr:  -21.9 


Figure  8.  Similar  to  Figure  7  but  derived  from  the  reciprocal  arra>  MAJO-South  The  prominent 
secondary  arrival  in  the  P  +  0  sec  window  has  a  phase  velocity  of  about  16  km/sec  and  azimuth 
of  about  260°. 


(a)  P  4  1,0  Hoc,  1  25  *  0.75  H/  (b)  I’  *  !  0  sue.  1  25  •  0.75  Hz 

Az:  19.1  Vol  800  Pwr:  16  7  A/  31  7  Vel:  9  99  Pwr:  218 


(c)  P  +  1.0  sec,  1.00  ±  0.50  Hz  (d)  P  +  7.0  sec.  1  00  +  0.50  Hz 

Az:  19.2  Vel:  8.89  Pwr  -15.8  Az:319  Vel:  9  97  Pwr:  20.9 


(e)  P  +  1  0  sec.  0.50  ±  0.30  Hz  (f)  P  +  7.0  sec.  0.50  ±  0  30  Hz 
Az:  1 9.8  Vel:  9  13  Pwr:  13.7  Az:  32  4  Vel:  9.48  Pwr:  18  1 


f  igure  9.  Similar  to  i;igurc  7  but  derived  from  lie  reciprocal  array  NA0.  The  prominent  secon¬ 
dary  arrival  in  the  P  +  7  sec  window  has  a  phase  velocity  of  about  1.8  km/sec  and  a/.imuth  of 
about  270°. 
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(.1)  P  <  10  r.nc,  1  25  *  0  75  It/  (b)  I1  i  70  sue.  1.25  ♦  0.75  Hz 

A/:  19.1  Vol.  8  90  Pwr:  16  7  A?  317  Vel:  9  99  Pwr:  21  8 


(c)  P  4-  1 .0  see,  1 .00  ±  0.50  Hz  (d)  P  +  7.0  sec,  1 .00  ±  0.50  Hz 


Az:  19  2  Vel:  8.89  Pwr:  15.8  Az:31.9  Vel  9  97  Pwi  20  9 


(e)  P  +  1 .0  sec.  0.50  ±  0.30  Hz  (f)  P  +  7.0  sec.  0.50  ±  0.30  Hz 

Az:19.8  Vel:  9  13  Pwr: -13.7  Az:  32.4  Vel:  9.48  Pwr: -18  1 


- 1  0  I  80  0  1  0  1  0  I  80  0  10 


Figure  9.  Similar  to  Figure  7  but  derived  from  the  reciprocal  array  NAO.  The  prominent  secon¬ 
dary  arrival  in  the  P  +  7  sec  window  has  a  phase  velocity  of  about  1.8  kni/scc  and  azimuth  of 
about  270°. 
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DISCUSSION 


F-k  analysis  of  NORESS  data  from  three  seismic  sources  provide  remarkably  consistent 
results  regarding  the  azimuthal  direction  and  phase  velocity  of  a  secondary  source  lying 
southwest  of  the  array,  in  the  region  of  Lake  Mjosa.  The  source-receiver  azimuthal  directions 
were  such  that  the  secondary  arrival  represents  back-scattering  for  the  Kazakh  explosion  and 
almost  normal  direction  of  scattering  for  the  other  two  events.  It  seems  therefore  that  the 
scattering  observed  in  this  study  is  not  strongly  dependent  on  the  direction  of  the  incident 
energy. 

Considering  the  local  structure  of  Yucca  Rat  as  a  basin  containing  low  velocity  tuffs  and 
alluvium,  Stead  and  Helmbergcr  (1988)  obtained  synthetics  that  matched  well  with  the 
observed  waveforms  and  also  explained  the  variation  of  coda  levels  with  shot  location  found 
by  Lynnes  and  Lay  (1988).  Gradual  basin  terminations  were  found  to  cause  the  greatest 
conversion  of  surface  wave  energy  into  tcleseismic  P  arrivals,  appearing  a  few  seconds  after 
P  and  pP.  Our  results  from  the  four  reciprocal  arrays  based  on  the  use  of  Yucca  Flat  explo¬ 
sions  also  indicate  secondary  sources  of  energy,  arriving  a  few  seconds  after  the  direct  P . 
Considering  the  seismic  velocities  in  this  region  (Hays  and  Murphy,  1971),  the  observed  low 
phase  velocity  suggests  these  arrivals  arc  due  to  Rayleigh  waves.  The  azimuthal  directions  of 
the  secondary  sources  in  our  study  show  large  variations.  Figure  11  shows  Ferguson’s  (198!) 
model  for  Yucca  Flat  basin,  derived  from  borehole  data  and  gravity  modeling,  along  with  the 
Yucca  Fault.  The  contours  indicate  depth  to  the  Tertiary-Paleozoic  contact  in  meters.  Loca¬ 
tions  of  the  39  shots  are  superposed  on  this  map.  The  secondary  arrival  observed  from  the 
MAJO-North  reciprocal  array  (Figure  7)  is  therefore  likely  to  be  due  to  Rayleigh  — »  P 
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scattering  near  the  northwestern  basin  boundary  of  the  Yucca  Valley.  Similarly,  our  results 
from  the  other  three  reciprocal  arrays  (Figures  X,  9,  and  10)  probably  indicate  arrivals  due  to 
Rayleigh  — >  P  scattering  near  other  basin  boundaries  of  the  Yucca  Valley  ( e.g .,  Lynnes  and 
Lay,  1989). 
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CONCLUSIONS 


Deterministic  locations  of  both  near-receiver  and  near-source  scattering  are  investigated 
by  using  several  f-k  techniques  and  polarization  analyses.  F-k  spectral  analyses  of  NORESS 
array  data  from  both  explosion  and  earthquake  sources  are  used  to  identify  and  locate  sources 
of  local  scattering.  Residual  seismograms  and  f-k  power  difference  plots  have  been  found  to 
be  useful  in  isolating  the  secondary  source  from  the  effects  of  the  primary  source.  The  most 
prominent  and  consistent  scattered  arrivals  have  been  identified  as  low-velocity  Rg  due  to  the 
scattering  of  incident  P  in  the  region  of  Lake  Mjosa  with  large  topographical  relief.  Four 
reciprocal  arrays  formed  by  using  single  station  records  of  closely  spaced  Yucca  Flat  explo¬ 
sions  are  used  to  investigate  near-source  scattering.  Considering  the  geological  structure  of  the 
Yucca  Valley,  the  secondary  sources  of  seismic  energy  appear  to  be  due  to  Rg  -»  P  scatter¬ 
ing  near  its  basin  boundaries. 
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